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1 Introduction 

One of the most surprising consequences of quantum mechanics is the entanglement of two or more 
distant particles. The two-particle entangled state was mathematically formulated by Schrodinger 
[1]: Consider a pure state for a system composed of two spatially separated subsystems, 

p=|*)(*|; |®> = X)c(a,6)|a)|6) (1) 

a,b 

where | a) and | b) are two sets of orthogonal vectors for subsystem 1 and 2. If c{a,b) does not 
factor into a product of the form f(a) x g(b), so that the state does not factor into a product 
state for subsystem 1 and 2 (e.g. p ± p\ 0p 2 ), the state is defined by Schrodinger as an “entangled 
state 11 . 

The classic example of a two-particle entangled state was suggested by Einstein, Podolsky, and 
Rosen in their famous 1935 gedankcnexperiment [2]: 


1^) = ~ c o)|«)|&) 

a, 6 


( 2 ) 


where c 0 is a constant. What is surprising about the entangled state (2) is the following: the 
measured value of an observable for either single subsystem is undeterminate. However , if one 
of the subsystems is measured to be at a certain value for that observable (the measured mine is 
certainly an eigen value) the other one is 100% determinate. This point can be easily seen from the 
delta function in state (2). Based on this unusual quantum behavior, EPR defined their “physical 
reality” and then asked the question: “Can Quantum- Mechanical Description of Physical Reality 
Be Considered Complete [3]?” One may not appreciate EPR’s criterian of physical reality and insist 
that “no elementary quantum phenomenon is a phenomenon until it is a recorded phenomenon 
[4], however, no one can ignore the unusual nonlocal behavior of state (2), especially considering 
when the measurements of subsystems 1 and 2 are space-like separated ( vents. 

Optical spontaneous parametric down conversion (SPDC) [5] [6] is the most effective mecha- 
nism to generate an EPR type entangled two-photon state. In SPDC, an optical beam, called the 
pump, is incident on a birefringent crystal. The pump is intense enough so that nonlinear effects 
lead to the conversion of pump photons into pairs of photons, historically called signal and i cl 1 pi . 
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The two-photon state generated from the SPDC crystal may be calculated from the standard 
theory (first order perturbation theory) to be [6], 

I* > = E *(“>• + «>i - “p + k . - M4 Mk s )) a t (u;,(k,)) |0) (3) 

s,t 

where u and k represent the frequency and the wave vector for signal (s), idler (i), and pump (p). 
The two delta functions in state (3) are usually explicitly written as phase matching conditions: 

u s + u, — Up , k s + k; = kp (4) 

Technically, the SPDC is said to be type I or type II, depending on whether the signal and 
idler beams have parallel or orthogonal polarization. The SPDC conversion efficiency is typically 
on the order of 10 -9 to 10~ 11 , depending on the SPDC nonlinear material. The signal and idler 
intensities are extremely low, only single photon detection devices can register them. It is clear 
that state (3) is an EPR type two-photon entangled state. The quantum entanglement nature of 
SPDC has been demonstrated in EPR-Bohm experiments and Bell’s inequality measurements [7]. 
The following two experiments were recently performed in our laboratory, which are more closely 
related to the original 1935 EPR gedankenexperiment. 

The first experiment is a two-photon optical imaging type experiment [8], which has been 
named “ghost image” by the physics community. The signal and idler beams of SPDC are sent 
in different directions, so that the detection of the signal and idler photons can be performed 
by two distant photon counting detectors. An aperture object (mask) is placed in front of the 
signal photon detector and illuminated by the signal beam through a convex lens. Surprisingly, 
an image of this aperture is observed in the idler beam, by scanning the idler photon detector 
in the transverse plane of the idler beam, if we are sure that the idler photon detector “catches” 
the “twin brother” of the signal , which can be easily performed by a coincidence measurement. 
This effect is even more striking when we found that the object-lens-image relationship satisfies 
the Gaussian thin lens equation. 

The second experiment demonstrates two-photon “ghost” interference-diffraction [9]. The 
experimental set up is similar to the image experiment, except that rather than a lens and an 
aperture it is a Young’s double-slit (or a single-slit) inserted into the path of the signal beam. We 
could not find any interference (or diffraction) pattern behind the slit. Surprisingly, an interference 
(or diffraction) pattern is observed when scanning the detector in the idler beam, if we are sure 
that the idler photon detector “catches” the “twin brother” of the signal. 

2 Two-photon “ghost” image experiment 

The experimental set-up is shown in fig. 1. The 351. lnm line of an Argon Ion laser is used to 
pump a nonlinear BBO (^-BaB 2 0 4 ) crystal which is cut at a degenerate Type-II phase matching 
angle to produce pairs of orthogonally polarized signal (e-ray of the BBO) and idler (o-ray of the 
BBO) photons [5]. The pairs emerge from the crystal near collinearly, with u, = = u p j 2, 

where Uj ( j = s,i,p) is the frequency of the signal, idler, and pump, respectively. The pump is 
then separated from the down conversion beam by a UV grade fused silica dispersion prism and 
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FIG.l. A schematic set-up for the two-photon “ghost” image experiment 

the remaining signal and idler beams are sent in different directions by a polarization beam splitting 
Thompson prism. The signal beam passes through a convex lens with a 400mm focal length and 
illuminates a chosen aperture (mask). As an example, one could choose the letters “UMBC”. 
Behind the aperture is the detector package D\, which consists of a 25 mm focal length collection 
lens in whose focal spot is a 0.8mm diameter dry ice cooled avalanche photodiode. The idler 
beam is met by detector package D 2 , which consists of a 0.5mm diameter multi-mode fiber whose 
output is mated with another dry ice cooled avalanche photodiode. The input tip of the fiber is 
scanned in the transverse plane by two encoder drivers, and the output pulses of each detector, 
which are operating in the Geiger mode, are sent to a coincidence counting circuit with a 1.8ns 
acceptance window. Both detectors are preceded by 83nm bandwidth spectral filters centered at 
the degenerate wavelength, 702. 2nm. 

By recording the coincidence counts as a function of the fiber tips’ transverse plane coordinates, 
we see the image of the chosen aperture (for example “UMBC”), as is reported in fig. 2. It is 



FIG. 2. a) A reproduction of the actual aperture “UMBC” placed in the signal 
beam. Note that the size of the letters is on the order of standard text, b) The 
image of “UMBC”: coincidence counts as a function of the fiber tip’s transverse plane 
coordinates. The scanning step size is 0.25mm. The data shown is a “slice” at the half 
maximum value, with no image enhancement. 
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interesting to note that while the size of the “UMBC” aperture inserted in the signal beam is 
only about 3.5mm X 7mm, the observed image measures 7mm X 14mm. We have therefore 
managed linear magnification by a factor of 2. Despite the completely different physical situation, 
the remarkable feature here is that the relationship between the focal length of the lens /, the 
aperture’s optical distance from the lens S, and the image’s optical distance from the lens (lens 
back through beamsplitter to BBO then along the idler beam to the image ) S' satisfies the Gaussian 
thin lens equation: 

111 /KX 

S + S' ~ f ^ 


In this experiment, we chose S = 600mm, and the twice magnified clear image was found when 
the fiber tip was in the plane with S' = 1200mm (see fig. 3). 

To understand this unusual phenomenon, we examine the quantum nature of the two-photon 
state produced in SPDC, which is entangled in momentum. The spatial distribution of the photon 
pairs is the result of the transverse components of the wave vector condition in equation (4) and 
Snell’s law upon exiting the crystal: 


k s sin ct s = ki sin cu — * u> s sin (3 S — sin /?, (6) 

where a s and a t are the scattering angles inside the crystal and 0 S and /?, are the exit angles 
of the signal and idler photons with respect to the k p direction. Therefore, near the degenerate 
case the photons constituting one pair are emitted at roughly equal, yet opposite, angles relative 
to the pump. Although the momentum of each photon is indeterminate, if one is measured at a 
certain value then the other is 100% determined. This then allows for a simple explanation of the 
experiment in terms of “usual” geometrical optics in the following manner: considering the action 
of the beamsplitter, we envision the crystal as a “hinge point” and “unfold” the schematic of fig. 1 
into that shown in fig. 3. Because of the equal angle requirement of equation (6), we see that 


fiber 
tip plane 



FIG. 3. A conceptual “unfolded” version of the schematic shown in fig. 1, which 
is helpful for understanding the physics. Although the placement of the lens and the 
detectors obey the Gaussian thin lens equation, it is important to remember that the 
geometric rays actually represent pairs of SPDC photons which propagate in different 
directions. 
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all photon pairs which result in a coincidence detection can be represented by straight lines (but 
keep in mind the different propagation directions) and therefore the image is well produced in 
coincidences when the aperture, lens, and fiber tip are located according to equation (5). In other 
words, the image is exactly the same as one would observe on a screen placed at the fiber tip if 
detector D x were replaced bv a point-like light source and the BBO crystal by a reflecting mirror 
[ 10 ]- 

3 “Ghost” interference-diffraction 

The schematic experimental set-up is illustrated in fig. 4. It is similar to the first experiment 
except that after the separation of signal and idler, the signal photon passes through a double-silt 
(or single-slit.) aperture and then travels about 1 m to meet a point-like photon counting detector 
I)\ (0.5mm in diameter). The idler photon travels a distance about 1.2m fiom Hb to the input 
tip of the optical fiber. In this experiment only the horizontal transverse coordinate, x 2 . of the 
fiber input, tip is scanned by an encoder driver. 


FIG. 4. A schematic set-up lor 
the two-photon “ghost” interference- 
diffraction experiment. The signal (e-ray 
of BBO) and idler (o ray of BBO) pho- 
ton pair is generated in nonlinear crys- 
tal BBO. The ultra violet pump beam 
is separated from the* down conversion 
beams by a I V grade lused silica disper- 
sion prism. BS is a beamsplitting Thomp- 
son prism for splitting the signal and idler 
beams to different directions. f\ and f 2 
are spectral filters with 702.2n.77i center 
wavelength and \0nm bandwidth. Both 
photon counting detectors D\ and D 2 are 
dry ice cooled avalanche photodiode op- 
erated in Geiger mode. 


Fig. 5 reports a typically observed double-slit interference-diffraction pattern. The coincidence 
counting rate is reported as a function of x 2 , which is obtained by scanning the detector D 2 (the 
fiber tip) in the idler beam, whereas the double-slit is in the signal beam. The Young’s double-silt 
has a slit-width a = 0.15mm and slit-distance d = 0.47mm. The interference period is measured 
to be 2.7 ± 0.2mm and the half-width of the envelope is estimated to be about 8mm. 
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Detector 2 position (mm) 


FIG. 5. Typically observed interference- 

diffraction pattern: the dependence of the co- 
incidences (per 400 sec) on the position of op- 
tical fiber tip of detector D 2 , which counts the 
idler photons, while the signal photons pass 
through a double-slit with a — 0.15mm and 
d = 0.47mm. The solid curve is calculated from 
equation (7). 



Detector 2 position (mm) 


FIG. 6. Two-photon diffraction: coincidence 
counts (per 400 sec) vs. the idler photon count- 
ing detector position. A single slit of width 
a = 0.4mm is in the signal beam. The solid 
curve is a fit from the theoretical calculations. 


By curve fittings, we conclude that the observation is a standard Young’s interference-diffraction 
pattern: 

R C { X ,) <x sine 2 (^) cos 2 . (7) 

The remarkable feature here is that 2 2 I s the distance from the slits plane, which is in the signal 
beam, back through BS to BBO crystal and then along the idler beam to the scanning fiber tip of 
detector Z) 2 (see fig. 4). The calculated interference period and the half-width of the sine function 
from equation (7) are 2.67mm and 8.4mm, respectively. Even though the interference-diffraction 
pattern is observed in coincidences, the single detector counting rates are both observed to be 
constant when scanning detector D\ or Z) 2 . Of course it is reasonable not to have any first order 
interference-diffraction in the single counting rate of Z) 2 , which is located in the “empty” idler 
beam. Of interest, however, is that the absence of the first order interference-diffraction structure 
in the single counting rate of D\, which is behind the double-slit, is mainly due to the divergence 
of the SPDC beam (> A /d). In other words, the “blurring out” of the first order interference 
fringes is due to the considerably large momentum uncertainty of a single SPDC photon. 

Furthermore, if Di is moved to an unsymmetrical point, which results in unequal distances to 
the two slits, the interference-diffraction pattern is observed to be simply shifted from the current 
symmetrical position to one side of a: 2 . This is quite mind boggling: imagine that there were a first 
order interference pattern behind the double-slit and D\ were moved to a completely destructive 
interference point (i.e. zero intensity at that point) and fixed there. Can we still observe the 
same interference pattern in the coincidences (same period, shape, and counting rate), except for 
a phase shift? 
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Fig. 6 reports a typical single-slit diffraction pattern. The slit-width a = 0.4mm. The pattern 
fits to the standard diffraction sine function, i.e., the “envelope” of equation (7), within reasonable 
experimental error. Here again z 2 is the unusual distance described in the above paragraphs. 

To explain this unusual phenomenon, we again present a simple quantum model where, sim- 
ilar to the “ghost” image experiment, the momentum of either single photon is undeterminate. 
However, if one is measured at a certain value the other one is determinate with unit probability. 
This important peculiarity selects the only possible optical paths in fig. 7, when one photon passes 
through the double-slit aperture while the other gets to D 2 . In the near degenerate case, we can 
simply treat the crystal as a reflecting mirror. 



FIG. 7: a) Simplified experimental scheme, and b) its “unfolded” version . 1 he 

overall optical path lengths between D x and D 2 along the upper (r^) and lower (r B ) 
paths, appearing in equation (11), are defined as: r A = r M + r A2 = r c i + r C2 , and 
r B = r B \ + r B 2 = rpi + rg 2 , where ra and roi are the optical path lengths from 
the slits C and D to the ith detector. 

The coincidence counting rate R c is determined by the probability P x2 of detecting a pair of 
photons by detectors D x and D 2 simultaneously. For SPDC, Pi 2 is proportional to the square of 
the second order correlation function (E { 2 +) E[ +) ) of the fields at points D x and D 2 (it thus plays 
the role of the two-photon effective wavefunction): 

Pi 2 = (E^E^E^E^) = |(4 +i p! +> )| 2 . (8) 

In equation (8) (...) = (<F| . . . |\I>), and I'P) is the four-mode state-vector of the SPDC field: 


|tf) = \vac) + e a} ate** + b\b]e i * B |«<*c) 


( 9 ) 


where t <C 1 is proportional to the pump field (classical) and the nonlinearity of the crystal, 4> A 
and <f> B are the phases of the pump field at A and B, and a+ (fc+) are the photon creation operators 
for the upper (lower) mode in fig. 7b (j = s, i). In terms of the Copenhagen interpretation one 
may say that the interference is due to the uncertainty in the birth-place (A or B in fig. 7) of a 
photon pair. 

In equation (8) the fields at the detectors are given by 


E[ +) = a s e ' kr + b s e ikTB ' ; 


4 +) = a t e' kr/ 


+ bit' 


kr F 


( 10 ) 
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where r Al (r Sl ) are the optical path lengths from region A (B) along the upper (lower) path to 
the ith detector. Substituting equations (9) and (10) into equation (8), 


Rc OC P 12 = c 2 \ e i(kr A +4>A) + ^(kr B +<t > B )]2 a 1 + cog _ rfi )] (1 1 ) 

where we assume <f> A = (j) B in the second line of equation (11) (although this is not a necessary 
condition to see the interference pattern, the transverse coherence of the pump beam at A and B is 
crucial). In equation (11) we defined the overall optical lengths between the detectors D\ and D 2 
along the upper and lower paths (see fig. 7b): r A = r M +r A2 = r cx +r C2 .r B = r H \+r B2 = r Dl +r D2 , 
where r Cl and rp; are the path lengths from the slits C and D to the ith detector. 

If the optical paths from the fixed detector D x to the two slits are equal, i.e., rc\ — r/p , and 
if z 2 d 2 j A, then r A — r B — rc 2 — r B2 = x 2 d/z 2 and equation (11) can be written as 

R c {x 2 ) oc cos 2 • (12) 

Equation (12) has the form of standard Young’s double-slit interference pattern. Here again 
z 2 is the unusual distance from the slits plane, which is in the signal beam , back through BS to the 
crystal and then along the idler beam to the scanning fiber tip of detector D 2 . 

If the optical paths from the fixed detector D\ to the two slits are unequal, i.e., rc 1 ^ t du 
the interference pattern will be shifted from the symmetrical form of equation (12) according 
to equation (11). This interesting phenomenon has been observed and reported following the 
discussion of fig. 5. 

To calculate the “ghost” diffraction effect of a single-slit such as shown in fig. 6, we need an 
integral of the effective two-photon wavefunction over the slit width: 

(l3 > 

where r(a- 0 , x 2 ) is the distance between points x 0 and x 2 ,x 0 belongs to the slit’s plane, and the 
inequality z 2 a 2 / A is assumed. 

Repeating the above calculations, the combined interference-diffraction coincidence counting 
rate for the double-slit case is given by 


R c {x 2 ) ex 


[ 2 dx 0 c~' kr(xo ' X2) 
J =r 


R c (x 2) oc sine 2 


/ x 2 r a\ 
V~A 7^) 


cos 2 



(M) 


which is exactly the same as equation (7) obtained from experimental data fittings. If the finite 
size of the detectors and the divergence of the pump are also taken into account by a convolution, 
the interference visibility will be reduced. These factors have been taken into account in the 
theoretical plots in figs. 5 and 6. 
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